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The  theoretical  analysis  required  to  determine  the 
electron  density,  electron  temperature,  excited  state  popu¬ 
lations  and  radiation  emission  of  monatomic  plasmas  is  des¬ 
cribed.  Available  information  on  inelastic  cross-sections 
and  radiative  emission  and  absorption  is  reviewed.  The  range 
of  physical  parameters  over  which  various  non-equilibrium 
effects  are  significant  is  examined.  Examples  of  the  appli¬ 
cation  of  this  theory  to  the  non- equilibrium-  magnetohydrtidynamic~ 
generator,  the  crossed-field  accelerator  and  the  shock  tube 
are  presented. 
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Introduction 


In  magnetohydrodynamics  it  is'  often  found  that  equilibrium 
c 

gas  properties  are  inadequate  in  describing  the  true  behavior  of  the  gas. 
Unfortunately  no  simple  prescription  has  been  evolved  which,  when  applied 
to  a  given  set  of  experimental  conditions,  determines  unambiguously  whether 
or  not  the  gas  is  in  equilibrium.  In  this  paper  a  summary  of  current  informa¬ 
tion  on  the  non-equilibrium  properties  of  partially  ionized  monatomic  gases 
is  given  from  which  a  rather  complex  prescription  can  be  formulated. 

If  one  or  more  of  the  following  criteria  are  met  the  gas  can 
be  said  to  be  out  of  equilibrium: 

1.  The  electron-ion-neutral  reaction  time  is  of  the  same >order 
as  or  longer  than  any  other  time  characteristic  of  the 
problem. 

2.  The  net  radiation  decay  rate  of  any  of  the  excited  states  is 
of  the  same  order  as  or  greater  than  the  net  collisional 
rate  of  populating  that  excited  state. 

3.  The  electron  temperature  is  significantly  different  from  the 
atom  and  ion  temperature. 

The  discussion  of  these  restrictions  is  limited  to  monatomic  gases  because 
of  the  absence  of  the  additional  complications  found  in  molecular  gases 
(rotation,  vibration,  dissociation,,  etc.)  and  also  because  only  for  monatomic 
gases  is  there  sufficient  information  on  collision  cross-sections  at  the 
present  time  to  make  use  of  these  criteria. 

It  has  been  established  previously^  that  electronic  excitation 
and  de-excitation  occur  primarily  by  electron  impact  with  atoms  and  by 
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radiative  emission  and  absorption.  Other  processes  such  as  dissociative 
recombination  and  atom  impact  excitation  do  occur,  but  they  are  insignificant 
in  monatomic  gases  when  the  electron  mole,  fraction  is  greater  than  10  "*  to 


A  semi-classical  method  for  computing  electron  impact  excitation 

and  de-excitation  cross-sections  has  been  developed  recently  by  Gryzinslci.'1' 

In  his  original  paper  Gryzinslci  shows  that  the  theory  agrees  well  with 

measurements  of  several  specific  excitations  and  ionization  cross-sections. 

7  8  9  ' 

The  theory  has  been  used  by  Bates,  et.al.  ’  and  by  Byron,  et.al.  to 
obtain  integrated  excitation  and  de-excitation  cross-sections  for  hydrogen¬ 
like  atoms.  Combining  these  cross-sections  with  available  radiative,  de- 

10  23 

excitation  probabilities  and  with  the  electron  energy  equation,  an 

electron-ion  recombination  theory  has  been  formulated  for  hydrogen  and  ■ 

9 

helium.  This  is  found  to  be  in  good  agreement  with  experimental  measure¬ 
ments  in  helium.*’’'1''1' 

In  the  following,  an  approximation  to  Gryzinslci 's  method  is  used 
to  compute  cross-sections  for  electronic  transitions  in  argon  and  potassium 
in  the  temperature  range  between  500°K  and  2000°K.  Coupling  these  with 
radiative  emission  and  absorption  probabilities,  the  rate  of  electron-ion 
recombination  is  determined  in  this  temperature  range.  The  rate  is  applied 
to  the  problem  of  producing  non-thermal  ionization  in  a  non-equilibrium 
mhd  generator.  Application  of  the  electron  energy  transfer  equation  to 
the  problem  of  producing  non-equilibrium  thermal  ionization  in  mhd  generators 
is  also  discussed.  The  influence  of  radiation  on  electron  temperature, 
population  of  the  excited  states,  and  electron  density  under  quasi-steady 
conditions  is  described.  Finally,  a  summary  is  given  of  the  application  of 
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the  reaction  'rate  equations  and  the  two-fluid  gasdynamic  flow  equations  to 
.compute  the  electron  density,  electron  temperature,  and  the  electrical 
conductivity  in  a  supersonic  nozzle  expansion  and  in  a  crossed  field  mhd 
accelerator . ^ 

Electronic  Transitions  by  Electron  Impact 

According  to  the  theory  of  Gryzinski ,  •*■  his  equation  (28)  is  used 
in  computing  the  cross-section  for  collisional  excitation  of  the  level  having 
binding  energy  E^,  to  that  having  binding  energy  ,  when  the  next 

higher  excited  state  has  a  binding  energy  -  U,.  The  cross-section  as 
a  function  of  electron  energy,.  E2,  is  then  integrated  over  a  Maxwell 
distribution  of  electron  energies.  The  result  for  hydrogen  is  given  in 
analytic  form  in  Ref.  9  and  shown  graphically  here  in  Fig.  3  as  a  function 
of  excited  state  binding  energy,  E  .  The  result  for  excited  states  of 
helium  in  the  binding  energy  range  shown  in  Fig.  3  is  approximately  twice 
that  of  hydrogen. 

'For  more  complex  atoms  such  as  argon  and  potassium,  a  great  number 
of  cross-sections  must  be  computed  and  summed  to  obtain  exact  results.  Noting 
that  the  cross-section  is  inversely  proportional  to  the  square  of  the 
energy  gap,  U,  between  the  initial  and  final  excited  states  (Eq.  23  and  26' 
of  Ref.  1)  an  approximation  is  immediately  suggested.  From  the  table  of 
atomic  energy  levels,  the  larger  energy  gaps  may  be  selected.  It  is  found 
that  the  energy  gaps  adjacent  to  the  large  gap,  U,  are  considerably  smaller 
than  U.  This  justifies  using  the  approximation  that  the  energy  levels 
above  and  below  the  gap  are  continuous.  It  is  also  found  that  these  closely 
spaced  energy  levels  extend  over  an  energy  range  greater  than  U  both  above 
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and  below  the  gap.  Therefore  if  the  electron  thermal  energy  is  of  the  order 
of  or  less  than  the  gap  energy,  U,  Eq.  26,  Ref.  1  is  a  good  approximation 
to  Eq.  28,  Ref.  1  in  computing  as  a  function  of  electron  energy. 

The  factor  E^/U)  found  in  eq.  26  and  27  for  this  case 

(U  +  E^  <■  E2)  is  reproduced  here  in  Fig.  1,  where  E^  is  the  kinetic  energy 

of  the  bound  electron  (the  same  as  the  binding  energy  for  a  single  excited 
electron),  U  the  energy  gap,  and  E2  the  kinetic  energy  of  the  impacting 
electron.  It  is  seen  that  a  linear  approximation  to  g^  vs  E2/U,  which 
greatly  simplifies  the  subsequent  integration  over  a  Maxwell  distribution 
of  electron  energies,  reproduces  g^  within  30%  over  the  range  1  +  .005  E^/U  < 
E2/U  £  2  .  1  E^/U.  If  we  make  the  assumption  that  kTg  lies  in  the 

range  .01  E^  $  kTg  ^  U  +  .1  E^,  the  linear  approximation  is  valid.  For 
example,  if  the  binding  energy  is  1  ev  and  the  gap  energy  is  0.2  ev,  the 
linear  approximation  is  valid  for  electron  temperatures  between  116°K  and 
3480°K.  The  integrated  cross-section  for  excitation  of  level  E^  to  E^  -  U 
or  higher  is  then  given  by 


o  /  oL> 

Q(u) c  f  J  ic'1 


JTkT* 

)  rwie 


.  (1) 


where  C~  is  the  collision  cross  section  (cm  ) ,  c  the  electron  speed 
(cm/sec),  Q(U)  is  obtained  from  eq.  (26),  Ref.  1,  f  is  the  Maxwell  distri- 

-1<6.  2  2 

bution  function,  CT  is  6.56  x  10  '  cm“  -  (electron  volt)  ,  B  is  the  slope 

of  the  linear  approximation  to  g^ >  U  is  the  gap  energy  (electron  volts),  k 
is  Boltzmann's  constant,  and  Tg  and  m^  are  the  electron  temperature  and 
mass.  The  slope,  B,  as  a  function  of  E^/U  is  reproduced  in  Fig.  2. 
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The  cross -sect! iii  given  in  eq.  1  is  for  a  transition  from  a 
particular  state  having  binding  energy  E^  to  any  state  having  binding  energy 
less  thin  or  equal  to  -  U.  The  total  rate  of  crossing  the  energy  gap 
U  in  the  upward  direction  is  actually  given  by  the  sum  of  the  contributions 
from  Ej.  and  all  levels  below  E^.  Using  the  approximation  that  the  levels 
below  E,  are  continuous,  the  total  rate  of  crossing  U  is 


where  and  are  the  number  density  and  degeneracy  of  the  states  having 
binding  energy  E^,  and  g  is  the  average  degeneracy  per  unit  energy  interval 
over  the  interval  U  below  E^. 

Using  the  principle  of  detailed  balancing  and  the  approximation 
of  continuous  levels  in  the  range  kTe  above  that  having  a  binding  energy, 

E-'  =  E^  -  U,  the  de-excitation  cross-section  from  all  levels  above  E“  and 
the  rate  of  crossing  the  gap  U  in  the  downward  direction  are  given  by 
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\Jhere  N?  and  g9  are  the  number  density  and  degeneracy  of  atoms  having  binding  energy 
E*  and  gU  is  the  total  degeneracy  of  states  having  binding  energy  between  E*  and  E*  - 
U.  At  equilibrium  N  is  related  to  the  free  electron  density  by 

CT 
D2 


N2  g  S.  '2  firm  kT 
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e  i 


,2  3/2  +E*/kT 

h  .  e  2 

•)  e.  N 
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where  h  is  planck1 s  constant  and  gg  and  are  the  degeneracy  of  the  free  electron 
and  the  ion. 

The  product  q-’c  has  been  determined  for  argon  and  potassium,  using  values 

13 

.for  U,  E.  ,  and  E*  derived  from  NBS  tables..  The  results  are  shown  graphically  in 


Fig.  3.  Note  that  these  cross-sections  are  unusually  lapge  ranging  between  10 

.  lf,-10  2 

and  10  cm  . 
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The  equilibrium  collisional  de-excitation  rate,  eqs.  (4)  and  (5),  normalized 
by  the  cube  of  the  electron  density,  is  shown  graphically  in  Fig.  4  for  hydrogen,  potas¬ 
sium,  and  argon  at  three  electron  temperatures.  Note  that  there  is  a  minimum  equilibrium 

9 

de-excitation  rate  for  each  gas  at  each  electron  temperature.  As  described  previously 
the  equilibrium  rate  of  de-excitation  at  the  minimum  point  determines  the  net  rate  of 
electron-ion  recombination,  the  ratio  of  the  latter  to  the  former  being  a  factor,  Y  , 
which  varies  between  1  and  1/4  (see  Ref.  9).  For  this  temperature  range,  Y-  1/3. 

_3 

The  resulting  electron-ion  recombination  rate  coefficient  N^  dN^/dt,  for 
argon,  potassium,  and  hydrogen  is  shown  graphically  as  a  function  of  electron  temperature 
in  Fig.  5.  Over  this  temperature  range  the  rate  for  helium  is  approximately  four  times 
that  for  hydrogen.  It  is  seen  that,  the  rate  for  potassium,  a  hydrogen-like  atom,  falls 
within  a  factor  of  2  of  that  for  hydrogen,  whereas  the  rate  for  argon,  an  atom  having 
very  closely  spaced  excited  states,  is  an  order  of  magnitude  larger  than  that  for 
hydrogen  over  quite  a  wide  range  of  temperature.  Note  also  that  the  recombination  rate 
changes  very  steeply  with  electron  temperature  in  this  range. 


Radiation 


At  low  electron  densities  radiative  transitions  contribute 
significantly  to  the  recombination  rate.  Also,  it  has  been  shown  that 
under  steady  state  conditions  radiation  shifts  the  distribution  of 
excited  states  away  from  the  equilibrium  distribution  and  may  lower 
the  electron  density  considerably  below  the  Saha  equilibrium  value. ^ 

The  spontaneous  emission  probabilities  for  hydrogen  are  found 
in  Refs.  10  and  14.  Values  for  helium  are  very  close  to  those  for 
hydrogen  except  for  low  lying  excited  states,  which  are  given  in  Ref.  15. 

Few  values  for  other  gases  such  as  potassium  and  argon  are  available. 

The  inclusion  of  radiation  in  the  ionization  and  recombination 

7 

theory  has  been  done  for  optically  thin  hydrogen  exactly  by  Bates,  et.al. 

9  ■ 

and  approximately  by  Byron,  et.al.  It  has  also  been  done  for  optically 

g 

thick  hydrogen  by  Ba'tes,  et.al.  The  resulting  recombination  rate  co- 

9 

efficients  for  the  optically  thin  case  are  shown  in  Fig.  5  for  hydrogen 
with  Ng  =  10^  and  10^  cm  ^ .  The  recombination  rate  coefficients  given 
here  for  hydrogen  are  as  much  as  a  factor  of  two  smaller  than  those  given 

by  Bates, ^  probably  because  of  the  simplifying  approximations  that  are 

,  9 
used . 

In  order  to  make  physically  realistic  computations  which  include 
radiative  effects  it  is  necessary  to  include  the  effect  of  re-absorption. 

This  can  be  done  in  a  general  way  by  considering  three  distinct  possibilities 
(1)  the  gas  is  optically  thick  only  for  transitions  from  excited  states 
to  the  ground  state,  (2)  it  is  optically  thick  for  all  line . radiation,  and 
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(3)  it  is  optically  thick  for  all  line  and  continuum  radiation.  Which 

of  these  models  if  any  applies  to  a  particular  physical  situation  depends 

on  each  of  the  absorption  coefficients, 
s 

2 

The  absorption  coefficient,  o/;  (cm  ),  of  an  atomic  line  is 
determined  by  the  oscillator  strength,  f,  of  the  line  and  the  line  shape. 
Only  the  strong  transitions  (f  -  1)  need  to  be  considered  since  these 
dominate  the  effect  of  radiation  on  excited  state  populations.  The  line 
shape  is  discussed  extensively  in  the  astrophysics  literature^  and  is 
thoroughly  reviewed  in  a  recent  article  by  Baranger . 17  The  principal 
sources  of  broadening  of  atomic  lines  in  partially  ionized  gases  of 
interest  here  are  Doppler  broadening  and  pressure  broadening. 

16 

Doppler  broadening  alone  can  be  treated  rather  simply, 
yielding  the  result  that  the  line  width,  aV  q  (sec"'1'),  and  the  absorp¬ 
tion  coefficient  at  line  center  are  given  by 

Av’  -  iTakr/M1  vjc  •  <YC  --  a  v** 

8  ,  (6) 

where  0  is  the  line  frequency,  M  the  atom  mass,  c  the  speed  of  light, 
m  the  electron  mass,  and  other  symbols  have  been  defined  previously.  The 
absorption  at  line  center  is  significant  here  because  emission  at  line 
center  contributes  most  to  changes  in  excited  state  populations.  The 

value  of  AVq/  for  plasma  temperatures  between  1000°IC  and  20,000°K 

-5  -6 

ranges  from  3  x  10  for  low  molecular  weight  gases  to  10  for  heavy 

gases  and  low  temperatures.  •  Using  f  1,  cYq  ranges  from  500/ V1'  to 

4  .  , 

1.7  x  10  /  Y Q-  Now  the  optical  depth  may  be  computed'if  the  number  density 
of  absorbing  atoms  is  known.  Using  an  optical  depth  of  1  cm  as  the  division 
between  optically  thick  and  optically  thin  problems,  the  division  occurs 
when  the- density  of  the  lower  state  is  in  the  range  h*  q/500  to  V*  /1..7  x  10^, 
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or  5  x  10  to  1.5  x  10  cm  for  10  ev  photons.  Since  the  ground 
state  population  in  magnetohydrodynamic  devices  is  usually  greater  than 
these  limits,  the  gas  is  optically  thick  for  transitions  from  an  upper 
state1  to  the  ground  state. 

The  excited  state  populations  (eq.  5)  are  of  the  order  of 

**  18  2 

10  and  the  emission  frequencies  for  which  the  lower  level  is  an 

excited  state  fall  in  the  range  1  to  4  ev.  Therefore  in  general  the  gas 

is  optically  thin  for  transitions  from  one  excited  state  to  another  if 

14  -3 

Ne  ^  10  cm  ,  whereas  it  is  optically  thick  for  these  transitions  if 

15  -3 

Ng  -5-10  cm  More  definitive  limits  than  these  can  be  set  only  by 
examining  each  physical  situation  in  detail. 

Pressure  broadening  exceeds  Doppler  broadening  of  low  excited 

15  -3 

states  at  electron  densities  of  the  order  of  10  cm  ,  and  lower  for 
higher  states.^  Techniques  are  now  available^  ’  ‘*'°  which  can  be  used 
to  compute  the  broadening  of  hydrogen  and  non-hydrogen  atoms,  but  few  speci¬ 
fic  examples  have  been  carried  out.  Considerable  work  in  the  future 
is  required  before  there  will  be  sufficient  information  on  line  broadening 
and  absorption  to  allow  accurate  quantitative  computations  of  the  effect 
of  radiation  on  gas  properties. 

It  should  be  pointed  out  here  that  molecular  band  emission 
covers  a  much  broader  wavelength  reginn  than  does  atomic  line  emission. 
Because  of  this  the  absorption  coefficient  is  much  lower  for  molecular 
emission  and  the  lower  state  population  required  to  produce  an  optical 
depth  of  1  cm  is  of  the  order  of  10^  cm  ^ .  Thus  molecular  band  emission 
is  generally  optically  thin. 
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The  free-free  and  free-bound  continuum  radiation  from  argon 

has  been  treated  theoretically  by  Pomerantz"^  and  studied  experimentally 
20 

by  Olsen.  Prom  these  results  it  is  found  that  an  absorption  length 

of  1  cm  for  the  continuum  in  the  visible  range  is  obtained  when  the  electron 

18  _  3 

density  is  of  the  order  of  10  cm  .  Since  in  mhd  devices  the  electron 

18  “3 

density  is  usually  much  smaller  than  10  cm  ,  the  gas  is  transparent 
to  continuum  radiation. 

Returning  to  the  three  classes  of  absorption  described  above, 
it  can  be  said  in  general  that  .the  gas  in  mhd  devices  is  optically  thick 
for  transitions  from  excited  states  to  the  ground  state,  optically  thin 
for  continuum  radiation,  and  the  optical  depth  for  transitions  from  one 
excited  state  to  another  may  be  of  the  order  of  1  cm  and  must  be  computed 
for  the  particular  operating  conditions. 

Non-Equilibrium  MHD  Generator 

The  electron-ion  recombination  rates  shoxm  in  Fig.  5  can  be 

4 

used  to  investigate  the  concept  of  non-thermal  ionization  in  mhd  generators. 

Because  mhd  generators  operate  efficiently  at  pressures  of  the  order  of  1 

14  -3 

atmosphere  or  greater,  electron  densities  must  be  of  the  order  of  10  cm 
or  greater-  in  order  to  reach  high  electrical  conductivity.  At  these  electron 
densities  and  in  the  temperature  range  typical  of  mhd  generators  radiation 
can  be  neglected  in  computing  the  recombination  rate  (see  Fig.  5).  It  is 
important  to  note  that  the  collisional  recombination  reaction  behaves 
like  a  three-body  process,  the  rate  being  proportional  to  the  cube  of  the 
electron  density,  not  to  the  square  of  the  electron  density  as  is  often 
assumed. 
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The  electrical  power  density  output  of  a  segmented  electrode 


generator  is 


r 


ar  cr  U 1  B 2 


where  the  load  factor  is  taken  as  1/2,  cr-  is  the  scalar  gas  conductivity, 

U  the  flow  velocity,  and  B  the  applied  magnetic  field  strength.  In  designing 
a  generator  the  magnetic  field  strength  is  chosen  to  produce  a  particular 
operating  value  of  the  Hall  parameter,  Af  =  10,  ^  ,  such  that  B  =  A  w /t  e 
xtfhere  U)  is  the  electron  cyclotron  frequency,  and  T  is  the  mean  time 
between  collisions  of  electrons  with  heavy  particles.  A  is  given  by 
'1  =  (N  (Qc)  ■*"  where  N  and'Q  are  the  number  density  and  collision 

cl  cl 

cross-section  for  the  heavy  particle  that  dominates  momentum  transfer 
collisions  with  electrons,  and  c  is  the  electron  mean  speed.  Then 


pe^  •=  ,  36"  cr  tt  ('A  >>!  ^  '&) 


(7) 


The  energy  required  to  replace  the  electrons  lost  by  recombination 
is  of  the  order  of 


P 

/o$ 


E' 


(7) 


where  E^on  is  the  ionization  potential  for  the  predominant  ion.  If  it 

is  assumed  that  P,  ^  0.1  P  for  efficient  operation,  a  limitation 

loss  gen  1 

on  the  operating  region  is'  given  by  the  inequality 


4;,  U,  *  .o^crU 


'-(  A?  c*  Va  8  c/e) 


(8) 
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Evidently  a  high  flow  velocity  and  high  Hall  parameter  are  favorable  whereas 
a  high  atom  density  is  unfavorable  because  it  must  be  accompanied  by  a  high 
electron  density  if  electrical  conductivity  is  not  to  -be  sacrificed. 

As  an  example,  helium,  at  a  stagnation  temperature  of  1800°K, 
operating  at  M  =  0.8,  T  =  1482°K,  U  =  2300  m/sec,  =  5,  N  =  3  x  1018  cm' 

JL  Q. 

is  chosen.  The  electron  neutral  cross-section  is  5.7  x  10  cm^  for  helium." 

14  -3 

To  approach  Spitzer  conductivity  the  electron  density  must  be  3.8  x  10  cm 

At  the  operating  temperature  =  9  x  10  ^  cm8/sec  for  helium  (four 

times  that  for  hydrogen,  see  Fig.  5).  Using  these  values,  the  loss  rate 

on  the  left  of  eq.  8  exceeds  the  allowed  supply  rate  on  the  right  by  a  factor 

of  2000.  One  operating  point  which  meets  the  restriction  of  eq.  8  is  that 

at  an  electron  density  that  is  lower  by  a  factor  of  j  2000,  or  a  value  of 
12  -3 

8.4  x  10  cm  ,  but  this  yields  an  electrical  conductivity  of  only  0.049 
mh<Vcm.  This  particular  example  does  not  appear  attractive  for  practical 
mhd  power  generation.  Further  detailed  analysis  may  yield  more  attractive 
operating  conditions. 

The  difficulty  encountered  here  can  be  eliminated  by  operating 

at  lower  pressures  and  higher  stagnation  temperatures.  Some  improvement 

is  to  be  found  using  different  gases,  for  example,,  a  potassium  additive 

to  helium  which  lowers  E.  by  a  factor  of  5  and  by  a  factor  of  2. 

ion  J  T 

The  use  of  some  other  noble  gas  such  as  argon  with  its  high  recombination 
coefficient  is  disadvantageous. 

An  interesting  possibility  is  that  the  recombination  process 
heats  the  electrons  to  a  temperature  substantially  higher  than  that  of  .  ■ 
the  gas.  This  reduces  the  recombination  rate  to  the  point  at  which  the 
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recombination  energy  is  equal  to  the  energy  transferred  by  elastic  collisions 

to  heavy  particles..  The- rate  of  transfer  of  energy  from  electrons  to  ions 

and  neutrals  is  given  by 
s 

to  --  3Hr&-ra)~  ( *4 4-,  i-  h 

(9) 

For  the  conditions  described  above  (and  assuming  that  radiation  losses 
are  unimportant)  the  electron  temperature  would  be  increased  to  3900°K, 
reducing  the  energy  loss  rate  to  5  times  that  of  the  electrical  power  density 
output. 

Operating  in  this  manner  the  "non-thermal  ionization"  mhd 

3 

generator  is  distinguished  from  the  "non-equilibrium  ionization"  generator 
only  by  the  manner  in  which  the  electrons  are  heated.  Joule  heating  by 
the  induced  current  is  used  in  the  latter,  whereas  an  electron  beam  or 
some  other  external  source  may  be  used  in  the  former. 

Electron  heating  in  non-equilibrium  mhd  generators  has  been 
2  3  5  26 

examined  in  some  detail.  *  ’  *  it  should  be  noted,  however,  that  a 
mistaken  impression  concerning  the  limits  of  the  mean  loss  factor,  £  , 
may  be  inferred  from  Ref.  3.  Following  eq.  (13)  it  is  stated  that  "The 
maximum  value  of  is  5/3,  corresponding  to  £  =  1,  K  “  0,  'and 

Y  =  5/3",  giving  the  impression  that  the  minimum  value  of  S  is  1. 

However,  if  a  small  amount  of  a  high  molecular  weight  easily  ionized  gas 
is  used  with  a  low. molecular  weight  gas  having  a  small  electron-neutral 
cross-section,  ^  may  be  much  less  than  1. 

If  two  components  are  chosen  for  the  working  fluid,  the  mixture 
ratio  that  minimizes  $  is  easily  determined  by  setting  the  derivative  of 
eq.  (1)  of  Ref.  3  equal  to  zero,  yielding  the  results 
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(10) 


where  M  is  the  molecular  weight,  Q  the  cross-section,  X  the  mole  fraction, 
and  the  subscript  2  represents  the  heavy,  large  cross-section  component 
of  the  gas.  For  C^/Q^  =  2000  which  is  roughly  the  ratio  of  the  Coulomb 
cross-section  of  cesium  to  the  ■  electron-neutral  cross-section  of  helium, 
the  optimum  mass. ratio  is  45,  also  close  to  that  for  cesium  and  helium. 

The  optimum  mole  fraction  of  cesium  ion  is  2.2%,  and  the  loss  factor, 

£  ,  under  these  conditions  is  0.09.  Using  eq.  15,  Ref.  3,  the  ratio 
of  electron  temperature  to  stagnation  temperature  is  equal  to  4.1  for  a 
segmented  electrode  generator  operating  with  this  mixture  at  M  =  0.8,  K  =  0. 

to  r  =  2. 

e  e 

To  summarize,  this  section,  calculated  values  of  the  electron-ion 
recombination  rate  are  so  high  that  "non-thermal  ionization"  for  mhd 
power  generation  is  probably  not  practical.  On  the  other  hand,  elevation 
of  the  electron  temperature  either  by  Joule  heating  or  by  external  methods 
or  both  appears  quite  practical,  particularly  when  applied  to  mixtures  of 
light  neutrals  with  heavy  ions. 
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QUASI-STEADY  EFFECTS 


Of  the  non- equilibrium  criteria  given  in  the  introduction, 
an  example  of  the  first  was  discussed  in  the  preceding  section.  The  second 
and  third  are  discussed  in  this  section.  For  optically  thin  hydrogenic 
atoms  under  quasi-steady  conditions.,  it  has  been  shown ^  that  a  significant 
reduction  in  the  populations  of  the  excited  states  and  the  electron  density 
below  the  equilibrium  values  based  on  statistical  mechanics  using  the  local 
value  of  the  electron  temperature  may  be  caused  by  radiation  from  the  ex¬ 
cited  states.  This  reduction  is  significant  if  the  net  radiative  de-excitation 
probability  of  an  excited  state  is  of  the  order  of  or  larger  than  the  net 
collisional  probability  of  exciting  the  state.  Application  of  this  criterion 
in  practice  requires  knowledge  of  the  collisional  excitation  cross-sections, 
the  radiative  decay  probabilities,  and  the  absorption  coefficients.  To  see 
the  range  of  physical  parameters  for  which  this,  effect  is  important,  hydrogen¬ 
like  atoms  subject  only  to  Doppler  broadening  will  be  examined. 


For  optically  thin  hydrogen  the  collisional  de-excitation  probability 

equals  the  radiative  de-excitation  probability  of  an  excited  state  having 

9  ■ 

principal  quantum  number  n  when 


X /  ‘/.s' 

/« 


-a  y 

IX 10  (n-,)  // 


(11) 


From  this  we  can  compute  the  electron  density  below  which  the  population  of 
level  n  may  differ  appreciably  from  equilibrium.  The  collisional  de-excitation 
probability  of  the  level  n  is  used  here  because  at  equilibrium  it  is  the  same  . 
as  the  probability  of  collisional  excitation  of  the  level  n. 


From  the  preceding  discussion  of  radiation,  the  population  of 

the  lower  state  below  which  a  sample  of  gas  1  cm  thick  is  optically  thin 

is  ‘‘■iven  by  N»  =  1/ c C  .  The  absorption  coefficient,  p/  ,  is  given  in 
°  ■■  J  J)  o  0 
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terms  of  the  line  width,  A^,  which  is  a  function  of  the  emission  frequency; 


Only  the  emission  from  a  level  n  to  the  next  lower  level  will  be 


considered  for  hydrogen,  as  it  is  assumed  that  the  gas-  temperature  is  so 

low  that  if  radiation  to  the  next  lower  level  is  absorbed,  emission  to 

all  levels  is  absorbed.  Then  V  =  "R,.  f  (w-i)  —  ft 

O'  y  l 

where  is  the  Rydberg  constant.  The  limit,  N  .  ,  is  then  given  by 


hi,  c  / 
nli  i  ,r  M 


m  R, 

h  e 


a  r  -a 

X  [  ( »->)  -  n 


(12) 


At  equilibrium,  the  population  of  the  level  n  -  1,  ,  is  related  to 

the  electron  density  by  eq.  5.  Then  eq.  12  yields  an  upper  limit  to  Ne, 

above  which  radiation  from  level  n  is  self-absorbed.  Eq.  11  yields  a  lower 

limit  to  below  xjhich  radiation  causes  the  population  of  level  n  to  be 

out  of  equilibrium.  Then  for  each  gas  temperature  a  level  nc  can  be  found 

at  which  these  two  restrictions  to  N  give  the  same  result.  It  can  then 

be  seen  that  if  the  electron  density  is  above  the  value  given  by  eq.  11 

using  the  value  of  nc  computed  from  eqs .  11,  12,  and  5,  all  of  the  excited 

states  will  be  nearly  in  equilibrium.  For  levels  below  n^  the  radiation 

is  trapped.  For  levels  above  nc,  the  collisional  excitation  rate  is  greater 

than  the  optically  thin  radiative  de-excitation  rate.  These  limiting  values 

of  n  and  N  are  tabulated  below  as  a  function  of  T  . 
c  e  e  • 


Table  I:  Electron  densities  below  which  non-equilibrium  steady- 
state  populations  of  the  excited  states  n  and  below 
may  occur  for  hydrogenic  atoms 


Te(  K)  •  43,500 


3 

5600 


1830 


5 

850 


Ne(cm"3)  3.7  x  1016  1.8  x  1014  8  x  1012 


8  x  10 


11 
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The  effect  of  pressure  broadening,  which  has  not  been  included 
in  computing  Table  I,  will  be  to  lower  the  temperatures  at  which  the 
limiting  values  of  Ng  given  in  Table  1  occur.  It  can  be  seen  that  many 
experimental  situations  arise  which  lead  to  non--equi  librium  steady  state 
populations  for  hydrogen.  For  other  gases  the  Doppler  line  broadening 
is  less  and  the  energy  levels  are  more  closely  spaced,  both  of  which  lower 
the  limiting  values  of  Ng. 

In  order  to  compute  the  gas  properties,  such  as  the  electron 
density  and  the  radiation  emission,  under  conditions  that  lead  to  non¬ 
equilibrium  excited  state  populations,  it  is  necessary  to  solve  the  set 
of  equations  that  includes  the  rate  of  populating  and  depopulating  each  level 
up  to  a  level  which  is  in  equilibrium  with  the  free  electron  density,  thus 

terminating  the  set  of  equations  with  eq.  5.  This  is  closely  related  to 

5  26 

the  theoretical  model  used  by  Ben  Daniel  *  to  determine  the  effect  of 
Joule  heating  of  electrons  on  electrical  conductivity. 

This  brings  us  to  the  third  non-equilibrium  effect  mentioned  in 
the  introduction.  The  electron  temperature  may  differ  from  the  atom 
temperature,  even  after  the  gas  relaxes  to  a  quasi-steady  value,  if  processes 
in  the  gas  supply  or  absorb  thermal  energy  of  the  electrons  more  rapidly 
than  it  can  be  transferred  from  electrons  to  ions  and  neutrals.  The 
following  discussion  will  be  focused  on  the  effect  of  radiation  on  the 
electron  temperature,  in  the  absence  of  electric  and  magnetic  fields. 

This  will  serve  to  delineate  the  regions  in  which  the  effect  of  radiation 
must  be  included  in  determining  the  electron  temperature. 
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For  this  purpose  radiation  can  be  divided  conveniently  into 
free-free  plus  free-bound  continuum  radiation  and  atomic  line  radiation, 
each  treated  separately.  The  electrons,  are  cooled  because  the  depletion 
of  excited  states  and  free  electrons  caused  by  radiative  transitions  is 
replenished  by  inelastic  collisional  excitation  and  ionization  by  electron 
impact.  This  results  in  a  lowering  of  the  electron  temperature  to  a  value 
which  produces  a  balance  between  the  inelastic  energy  loss  rate  (equal 
to  the  radiation  loss  rate)  and  the  elastic  energy  transfer  rate  between 
electrons  and  heavy  particles. 


The  continuum  radiation  is  proportional  to  the  square  of  the 
19  22 

electron  density,  ’  as  is  the  rate  of  transfer  of  energy  from  electrons 
23 

to  ions.  Since"  the  former  is  independent  of  ion  mass  and  the  latter 

inversely  proportional  to  ion  mass,  the  largest  temperature  difference 

is  expected  for  the  heaviest  iions.  In  the  event  that  elastic  energy 

transfer  from  neutrals  to  electrons  is  dominant,  the  temperature  difference 

will  be  smaller  than  that  computed  on  the  assumption  that  ion-electron 

elastic  energy  transfer  is  dominant.  If  this  is  done  by  equating  the  total 
22 

continuum  radiation  to  eq.  11,  Ref.  23,  the  electron  temperature  is 
given  by 


4  ✓  h)  _ _ . 

(  ri  ~  1 '  "*■  i/T  JU  [  <?(kre)  "'/zrrtfe  e  *J 


where  aV  is  the  frequency  range  over  which  the  continuum  extends.  There 

20  22 

is  recent  experimental  evidence  that  the  Unsold  theory  underestimates 

the  continuum  of  argon  by  a  factor  of  5.  However,  the  total  continuum 

energy  loss  rate  is  within  a  factor  of  two  of  the  experimental  value  if 

«  1.5  x  10^  sec  ^  is  used.  Using  this  it  is  found  from  eq.  13  that 

T  /T  =  1.02  in  xenon,  showing  that  the  continuum  radiation  is  not  sufficient 
a  e 


-18- 


under  any  circumstances  to  cause  the  electron  temperature  to  differ 
appreciably  from  the  gas  temperature. 


Turning  now  to  the  effect  of. atomic  line  radiation  on  the  electron 

s 

temperature,  we  have  for  the  energy  radiated  per  unit  volume  per  second, 


•  '*•  %  Q 

w/i t-6  j  o 


where  y  is  the  emission  frequency,  A^  the  transition  probability,  and 

N.  the  density  of  the  state  j.  The  sum  over  j  is  limited  to  optically  thin 

emission  lines.-  For  strongly  allowed  transitions  the  emission  probability 

2  2  2  3 

is  proportional  to  the  square  of  the  frequency^,  A..  -  8  tr  e  /g^mc  . 

The  upper  level,  N,,  is  related  to  the  lower  level  by  N./N*  =  g.exp(-h  v>?  ,/kT  )/g 
,  .  J  J  •*  3  ’  3  e. 

If  the  gas  is  to  be  optically  thin  the  lower  level  is  limited  to  the  value 
given  by  eq.  12  with  £  (n  -  1)  2  -  n  2  J  replaced  by  h  Then 


Q.  -  I  I  ^  exp  / -  b  *  / kTe) 

c3  I  irH  J  /  J  (15) 

It  is  seen  that  Q  has  a  maximum  at.h  V.  *  V  kT  .  Substituting  this  value 
J  J  e 

into  eq.  15  and  retaining  only  one  level,  j,  as  representative  of  the  sum, 

23 

eq.  14  can  be  equated  to  the  elastic  energy  transfer  rate  to  give  for 
the  temperature  difference 


7k  _  /  -  5-1 J KkTa/r»J  l-KTi)  /h  c  e  He 

Te  r  , 
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Using  potassium  at  Tg  “  3000°K  as  a  numerical  example,  Ta/Tg  -  1.2  when 
13  -3 

N  =  1.6  x  10  cm  .  This  is  a  lower  limit  to  N  for  this  temperature 
6  © 

difference  because  only  one  line  was  included  in  the  sum  and  also  because 
pressure  broadening  may  appreciably  reduce  the  absorption  coefficient. 

A  lowering  of  the  electron  temperature  below,  the  gas  temperature 
under  quasi-steady  conditions  has  been  reported  in  the  shock  tube  measure- 

24 

ments  of  excitation  temperatures  of  chromium  done  by  Charatis  and  Wilkerson. 

The  measured  electron  temperature  was  of  the  order  of  6000  K  for  an  atom 
temperature  of  10,000°K.  However,  in  a  more  recent  article  Charatis  and 
Wilkerson  have  revised  the  radiative  transition  probabilities  used  previously 
by  computing  them  from  their  spectroscopic  measurements.  They  now  find  ( 

that  the  electron  temperature  is  within  3%  of  the  gas  temperature.  Considerable 
caution  must  be  exercised  in  using  one  experiment  to  show  that  the  gas  is 
in  equilibrium  and  also  to  obtain  radiative  transition  probabilities.  Further 
experimental  data  is  needed  in  order  to  determine  with  certainty  the  influence 
of  line  emission  on  the  electron  temperature.  Note  also  that  the  introduction 
of  stable  molecular  species  will  greatly  increase  the  radiation  energy  loss 
and  may  lower  the  electron  temperature  significantly. 

To  summarize  this  section^  the  range  of  physical  parameters  for 
which  non-equilibrium  populations  of  excited  states  may  occur  has  been 
determined-  and  generally  falls  outside  the  operating  range  of  mhd  power 
generators,  but  does  ,  .  include  the  operating  range  of  many  propulsion  de¬ 

vices.  The  continuum  radiation  has  been  shown  to  have  negligible  effect  on 
the  electron  temperature.  The  range  of  parameters  for  which  line  radiation 
significantly  lowers  the  electron  temperature  has  been  given  and  also  gener¬ 
ally  falls  outside  the  operating  range  of  mhd  power  generators,  but  includes 
the  operating  range  of  many  propulsion  devices. 


Nozzle  Expansion  and  Crossed  Field  Accelerator 


The  reaction  rate  theory  described  here  has  been  combined  with 

two  ^luid  one-dimensional  compressible  flow  theory  and  solved  on  an  IBM  7090 

to  determine  the  properties  of  argon  and  of  helium  in  the  deLaval  nozzle 

12 

expansion  of  the  thermal  arc  heated  gas.  It  has  also  been  used  to  compute’ 

electron  density,  electron  temperature,  electrical  conductivity,  gas 

.temperature,  and  enthalpy  in  a  supersonic  free  jet  crossed  field  accelerator. 

The  accelerator  performance  determined  theoretically  is  in  reasonable  agree- 

12 

ment  with  experimental  measurements.  A  few  typical  results  are  shown 
graphically  here  in  Figs.  6,  7,  and  8  and  are  discussed  briefly  below. 

For  details  see  Ref.  12. 

The  degree  of  ionization  of  helium  in  a  deLaval  nozzle  expansion 

is  shown  in  Fig.  6.  It  is  assumed  that  the  gas  is  in  equilibrium  initially 

at  the  nozzle  entrance.  Two  solutions  are  shown,  one  assuming  that  radiation 

is  absorbed,  the  other  assuming  that  all  radiation  is  emitted.  The  electron 

density  for  the  optically  thin  case  decreases  even  in  the  nozzle  entrance, 

tending  toward  the  quasi-steady  value  of  electron  density  for  the  case  of 

optically  thin  helium  at  the  stagnation  temperature.  Based  on  estimates 

of  the  absorption  given  above,  the  optically  thick  case  is  applicable  to 

most  experimental  conditions.  •  In  both  cases  it  is  seen  that  the  amount  of 

recombination  is  small.  It  is  also  evident  that  the  approximation  suggested 

27 

by  Bray  which  has  been  successfully  applied  to  atomic  recombination  in 
nozzle  expansions  is' not  valid  for  electron-ion  recombination.- 

Results  of  calculations  of  the  properties  of  argon  in  a  supersonic 
free  jet  crossed-field  accelerator  are  shown  in  Figs.  7  and  8.  The  starting 
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conditions,  noted  in  the  figure,  are  those  obtained  from  the  nozzle 

computation  described  above  and  applied  to  argon.  A  step  function  increase 

in  electric  and  magnetic  field  at  x  =  0  is  assumed.  In  this  calculation 

the  effect  of  radiation  on  the  reaction  rates  is  neglected.  This  is  not 

justified  for  hydrogenic  atoms,  but  is  for  argon  because  of  the  closely 

spaced  excited  electronic  levels.  It  has  also  been  shown  that  radiation 

energy  losses  are  negligible  in  the  crossed  field  accelerator  operating 

12 

under  the  conditions  of  the  computations  presented  here.  In  Fig.  7  it 
is  seen  that  the  electron  temperature  rises  sharply  to  about  16,000°K, 
followed  by  a  region  of  increasing  electron  density.  An  initial- "dip  in 
the  electrical  current  is  caused  by  an  increase  in  the  Hall  parameter,  due 
to  the  increase  in  electron  temperature  and  corresponding  decrease  in 
cross-section.  The  electrical  current  reaches  zero  at  the  point  at  which 
the  gas  has  been  accelerated  to  the  E/B  velocity.  Very  little  electron- 
ion  recombination  occurs  in  the  exhaust  of  the  accelerator  because  of  the 
low  electron  density  and  high  electron  temperature. 
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Conclusions 

The  basis  for  a  satisfactory  electrori-ion  reaction  rate  theory 

has  b,een  presented  and  used  to  indicate  the  range  of  physical  parameters 

for  which  significant  non-equilibrium  pehenomena  may  be  encountered. 

Illustrations  of  the  application  of  the  theory  have  been  given  from  which 
'  '  , 
it  appears  that  the  use  of  non-thermal  ionization,  that  is  ionization 

greater  than  the  Saha  equilibrium  value  at  the  electron  temperature,  in 

mhd  generators  is  not  practical,  and  that  electron  heating  either  by 

internal  current  or  by  external  means  in  non-equilibrium  mhd  generators 

is  most  attractive  when  applied  to  mixtures  of  heavy  ions  and  light  neutrals. 

It  was  shown  that  the  influence  of  line  radiation  on  the  electron  temperature 

is  significant  at  low  electron  densities  although  continuum  radiation  .is 

not.  The  electron  density  and  temperature  in  an  arc  jet  nozzle  expansion 

and  a  crossed  field  accelerator  are  found  to  be  far  from  equilibrium 

which  shows  that  a  complete  reaction  rate  theory  is  required  in  computing 

gas  properties  under  these  conditions. 
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‘Figure  Captions 


1. 


2. 


The  function  g.(E_/l';  E^/U)  used  in  Gryzinski's  semi-classical  theory 
of  Electronic  inelastic  cross-sections, 
s 


Slope,  B,  of  the  portion  of  that  is  a  linear  function  of  E2/U. 


3.  The  averaged  product  of  the  electron  speed  times  the.  collisional 
de-excititation  cross-section  of  the  excited  state  having  binding 
energy  E”  as  a  function  of  E*  for  several  gases. 

4.  De-excitation  rate  of  excited  states  normalized  by  the  cube  of  the 

electron  density  as  a  function  of  i:he  binding  energy.  These  rates  are 
computed  assuming  that  the  excited  states  are  maintained  in  equilibrium 
with  the  free  electron  density  by  inelastic  collisions.  v 

5.  Theoretical  electron-ion  recombination  rates  normalized  by  the  cube  of 
the  electron  density  for  several  gases. 


6.  Degree  of  ionization  of  helium  in  a  deLaval  nozzle  expansion  from  a  thermal 
arc. 

7.  Gas  properties  of  argon  in  a  supersonic  free  jet  crossed  field  accelerator 
using  a  constant  pressure  two  fluid  non-equilibrium  compressible  flow  theory. 

8.  Gas  properties,  of  argon  in  a  supersonic  free,  jet  crossed  field  accelerator 
using  a  constant  pressure  two  fluid  non-equilibrium  compressible  flow  theory. 
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